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seen  in  brine  pocket  structure  between  the  years.  Values  of  ef.  as  a  function  of  i&  were 
slightly  lower  on  the  average  at  9.50  GHz  than  at  4.80,  and  most  values  of  e^.  agreed 
with  the  results  of  others  when  their  temperature  and  salinity  values  were  converted 
to  values  of  Changes  in  £‘r  due  to  thermal  modification  (laboratory  warming)  seem 
more  strongly  related  to  a  gain  in  air  volume  from  brine  drainage  than  to  the  redistri¬ 
bution  of  brine  inclusions.  The  in-situ  measurements  showed  extremely  high  attenuation 
for  the  young  (<  12  cm)  brine-rich  ice.  Good  agreement  was  found  between  data  for 
the  more  desalinated  samples  and  the  theoretical  values  predicted  by  a  previously 
proposed  dielectric  mixing  model.  This  model  was  modified  by  accounting  for  the 
brine  pocket  geometry  observed  in  thin  sections,  and  also  by  including  a  bulk  conduc¬ 
tivity  term  to  account  for  the  observed  loss. 
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PREFACE 
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and  Engineering  Laboratory. 

This  work  was  performed  as  part  of  a  saline  ice  research  program  conducted  for  the  Office 
of  Naval  Research  (ONR)  under  Contract  No.  N000148WRM24001. 

The  assistance  of  Allan  Delaney  in  many  aspects  of  this  project  is  gratefully  appreciated. 
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Microwave  and  Structural  Properties  of  Saline  Ice 

ANTHONY  J.  GOW,  STEVEN  A.  ARCONE  AND  SETH  G.  MCGREW 

Part  I.  Structural  Properties 


INTRODUCTION 

Despite  widespread  imaging  of  sea  ice  in  the 
Arctic  and  Antarctic  by  various  remote  sensing 
techniques,  only  limited  attempts  have  been  made 
to  obtain  precise  correlation  of  the  electrical  prop¬ 
erties  of  sea  ice  with  its  structural  and  salinity 
characteristics.  Such  correlations  (or  calibrations) 
can  be  accomplished  either  by  conducting  ground 
truth  measurements  simultaneously  with  imaging 
overflights  or  by  performing  similar  scale  model 
tests  on  simulated  sea  ice  grown  in  a  test  tank  or 
pool.  To  this  end,  an  outdoor  pool  was  construct¬ 
ed  at  the  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory  (CRREL)  and  a  multidis¬ 
ciplinary  investigation  undertaken  to  examine  the 
effect  of  the  salinity  and  structural  characteristics 
of  sea  ice  on  reflectivity  and  emi"sivity  properties. 
The  purpose  of  Part  I  of  this  report  is  to  document 
the  salinity  and  structural  characteristics  of  the 
simulated  sea  ice  (saline  ice)  used  for  making 
measurements  of  dielectric  properties  under  care¬ 
fully  controlled  conditions  in  the  laboratory.  Re¬ 
sults  of  the  dielectric  measurements  are  reported 
by  Arcone  and  McGrew  in  Part  II  of  this  report — 
Microwave  Properties.  This  report  is  based  in  part 
on  results  presented  by  Arcone  and  McGrew  (1985) 
and  Gow  (1985)  at  the  1985  International  Geosci¬ 
ence  and  Remote  Sensing  Symposium  (IGARSS 
‘85)  held  at  the  University  of  Massachusetts, 
Amherst. 

EXPERIMENTAL  SET-UP 

Ice  sheets  were  grown  in  a  pool,  measuring 
12.2x5.2  m,  that  was  excavated  to  a  maximum 
depth  of  1.5  m  and  covered  with  a  thick  neoprene 
liner  (Fig.  1).  Concrete  pavers  were  emplaced 


around  the  edge  of  the  pool  to  serve  as  a  work 
pad.  The  rest  of  the  facility  consisted  of  a  large 
tent  that  could  be  rolled  on  and  off  the  pool  as 
weather  conditions  dictated,  and  a  refrigeration 
unit  that  was  used  primarily  to  offset  any  melting 
of  the  ice  during  periods  of  elevated  surface  air 
temperature.  A  space  heater  was  also  installed  to 
suppress  excessive  ice  growth  during  periods  of 
very  cold  weather.  Other  investigations  concurrent 
with  CRREL’s  efforts  (Grenfell  and  Comiso  1986, 
Swift  et  al.  1986)  utilized  remote  sensing  antennas 
and  other  instruments  mounted  on  a  steel  gantry 
that  spanned  the  pool  and  could  be  moved  across 
the  pavers  to  permit  imaging  measurements  on 
any  part  of  the  ice  sheet.  In  total,  personnel  from 
the  Universities  of  Kansas,  Washington,  Massa¬ 
chusetts  and  Wisconsin,  the  Environmental  Re¬ 
search  Institute  of  Michigan  (ER1M),  Massachu¬ 
setts  Institute  of  Technology,  NASA  and  the  U.S. 
Navy  laboratories,  NORDA  and  NUSC,  also  partic¬ 
ipated  in  the  experiments  that  included  radiom¬ 
eter,  radar  and  acoustical  backscatter  and  dielec¬ 
tric  property  measurements. 

The  pool  was  filled  with  water  that  was  raised  to 
the  desired  salinity  (23-25%o)  by  adding  granulat¬ 
ed  sea  salt  and  mixing  it  thoroughly  with  circulat¬ 
ing  pumps. 

Full-scale  testing  was  begun  on  19  December 
1983  when  the  pool  was  allowed  to  freeze  spontan¬ 
eously  without  seeding.  Growth  of  this  ice  sheet 
resulted  in  the  formation  of  massive  prismatic 
crystals  with  c-axes  that  reverted  slowly  to  hori¬ 
zontal  orientations  at  the  bottom.  At  this  scale 
such  ice  does  not  accurately  reflect  the  textures  of 
naturally  frozen  sea  water.  Accordingly,  this  ice 
sheet  was  broken  up  and  the  blocks  entirely  re¬ 
moved  in  readiness  for  a  second  ice  sheet.  On  17 
January  1984,  the  surface  of  the  pool  was  spray 


b.  Experimental  ice  sheet  with  tent  rolled  of  t 
Figure  I.  Saline  ice  test  facility  at  (HE 1 1 


seeded  in  an  attempt  to  stimulate  growth  of  more 
typically  columnar  ice  crystals.  This  technique, 
previously  employed  by  the  Gow  (1986)  to  pro¬ 
mote  columnar  crystal  growth  in  simulated  fresh¬ 
water  ice,  also  succeeded  for  the  saline  water  in  the 
pool.  The  resultant  ice  sheet  was  allowed  to  con¬ 


tinue  growing  for  36  davs  ot  testing  during  which 
the  structural  and  salimtv  characteristics  were 
monitored  at  regular  intervals.  The  crystalline  tex¬ 
ture  of  the  ice  at  all  stages  of  its  growth  was  found 
to  closely  scale  that  of  normal  congelation  ice. 
Congelation  ice  forms  bv  direct  freezing  of  sea 


water  to  the  underside  of  the  sheet;  such  a  process 
typically  yields  vertically  elongated  crystals  with 
their  c-axes  oriented  substantially  in  the  horizontal 
plane  (Weeks  and  Ackley  1982). 

Experiments  were  continued  on  an  expanded 
scale  during  the  1984-85  winter  when  three  ice 
sheets  were  grown  and  tested.  These  included  one 
ice  sheet  that  was  divided  into  seeded  and  unseed¬ 
ed  portions,  and  two  seeded  ice  sheets,  the  final 
one  of  which  we  succeeded  in  converting  into  a 
sheet  with  structural  and  salinity  properties  transi¬ 
tional  to  those  of  second-year  ice. 


ANALYTICAL  TECHNIQUES 

Once  freezing  had  begun,  each  ice  sheet  was 
monitored  closely  for  salinity  and  structure 
throughout  its  entire  growth  history.  To  limit  the 
effects  of  brine  drainage,  salinity  measurements 
were  usually  made  on  specimens  prepared  on  site 
as  soon  as  a  sample  was  sawed  or  cored  from  the 
ice  sheet.  Often  additional  measurements  were  run 
on  a  second  core  or  sample  to  determine  if  salinity 
differed  significantly  on  different  parts  of  an  ice 
sheet.  Results  of  these  measurements  generally 
agreed  within  0.3VM  at  corresponding  levels  in  the 
ice  sheet.  Salinities  were  measured  with  a  Beck¬ 
man  Solubridge  RB-5  that  was  frequently  cali¬ 
brated  against  solutions  prepared  from  Copen¬ 
hagen  Standard  Seawater  (chlorinity  =  19.373%o). 
The  precision  of  these  measurements  was  estimat¬ 
ed  at  ±0.2%o.  Data  were  either  plotted  conven¬ 
tionally  as  a  function  of  depth  or  as  cumulative 
depth  (depth  x  salinity)  values  vs  depth. 

Crystal  structure  studies  were  made  with  the  aid 
of  thin  sections  cut  vertically  and  horizontally 
from  blocks  sawn  from  the  ice  sheet,  thereby  fur¬ 
nishing  a  three-dimensional  view  of  model  sea  ice 
structure.  Sections  were  prepared  on  a  histological 
microtome  adapted  for  use  with  ice.  Details  of  this 
technique  as  applied  to  sea  ice  are  given  in  Gow 
and  Weeks  (1977).  Sections  were  sliced  to  a  thick¬ 
ness  of  0.5  mm  or  less  and  then  photographed  be¬ 
tween  crossed  Polaroids  at  two  magnifications:  at 
natural  (1:1)  scale  with  a  10-xl2.5-cm  press 
camera  and  at  7  x  magnification  with  a  12.5- x 
17.5-cm  bellows  extension  camera.  The  latter  for¬ 
mat  yields  large  photos  for  examining  the  finer  de¬ 
tails  of  ice  platelet  and  brine  inclusion  characteris¬ 
tics  in  individual  saline  ice  crystals. 


RESULTS  AND  DISCUSSION 

The  natural  freezing  of  seawater  generally  leads 
to  growth  of  vertically  elongated  columnar  crys¬ 
tals  in  sheets  that  reach  2  m  or  so  in  thickness. 
Closer  inspection  of  the  columnar  crystals  invaria¬ 
bly  reveals  a  substructure  of  long,  vertical  plates 
separated  by  parallel  layers  of  brine  inclusions  (see 
Fig.  2).  Such  vertical  sandwich  substructure  re¬ 
flects  the  process  by  which  excess  brine  at  the  ice/ 
water  interface  is  systematically  incorporated  into 
the  ice.  The  spacing  of  these  plates  (plate  width) 
can  vary  from  a  few  tenths  of  a  millimeter  to  1  mm 
and  depends  mainly  on  the  speed  of  growth;  the 
faster  the  freezing,  the  narrower  the  plate  spacing 
and  the  greater  the  salinity.  Subsequent  changes  in 
this  substructure  occur  mainly  in  response  to  tem¬ 
perature  changes  in  the  ice.  Brine  inclusions  are 
particularly  sensitive  in  this  regard.  Even  small 
changes  in  thermal  regime,  resulting  from  day-to- 
day  changes  in  surface  air  temperatures,  can  lead 
to  significant  changes  in  the  geometry  of  the  inclu¬ 
sions  and  the  concentration  of  the  entrapped  brine. 


1983-84  experiments 

Salinity  and  structure  characteristics  of  the 
seeded  ice  sheet  are  documented  in  Figures  3 
through  6.  Magnified  sections  of  crystal  substruc¬ 
ture  in  Figure  6  show  how  well  this  simulated  sea 
ice  matches  the  naturally  frozen  product  (Fig.  2). 
This  ice  sheet  had  attained  a  thickness  of  32.5  cm 
before  measurements  were  terminated  on  22  Feb¬ 
ruary  1984,  36  days  after  ice  growth  was  initiated. 
A  maximum  bulk  salinity  of  80/M  was  measured  at 
the  end  of  three  days  of  growth  (see  Fig.  3),  with 
salinities  then  falling  progressively  to  4%o  after  20 
days  of  growth  and  to  2%o  after  36  days.  This  de¬ 
salination,  particularly  that  for  the  interval  be¬ 
tween  20  and  36  days,  can  be  correlated  with  per¬ 
iods  of  elevated  air  temperature  that  led  to  flood¬ 
ing  and  refreezing  of  5  cm  of  freshwater  on  top  of 
the  ice  sheet  during  the  latter  stages  of  the  experi¬ 
ment.  These  changes  are  also  demonstrated  in  the 
cumulative  salinity  profiles  presented  in  Figure  3b. 

Figure  4  shows  the  overall  nature  of  the  in-situ 
ice  structure  observed  in  vertical  and  horizontal 
thin  sections  of  15-day-old  ice  photographed  be¬ 
tween  crossed  Polaroids.  Note  the  columnar  tex¬ 
ture  of  crystals  in  the  vertical  section  (some  of 
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Figure  2.  Horizontal  thin  section  of  arctic  sea  ice,  showing  crystal  substructure 
composed  of  ice  platelets  and  brine  layers.  The  crystallographic  c-axis  (arrow)  lies  in 
the  plane  of  the  section  and  is  oriented  normal  to  the  platy  structure  of  individual  crystals. 
Scale  subdivisions  are  I  mm. 


a.  Bulk  salinity. 


20  40 

Age  (days) 


6  Days  /  / 


<  /  / 
/  A 


(0  20 

Ice  Thickness  (cm) 

b.  Cumulative  salinitv. 


Figure  3.  Bulk  salinity  and  cumulative  salinity 
( salinity  x  ice  thickness  increment)  measured  as  a 
function  of  the  age  and  thickness  of  the  ice,  respec¬ 
tively.  The  ice  sheet  was  grown  during  January  and  Feb¬ 
ruary  1 984.  In  b  the  profile  for  an  ice  sheet  gro  wing  with  a 
uniform  salinity  of  6°/m  is  shown  for  comparison. 
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Figure  4.  Vertical  and  horizontal  thin  sections  of  structure  of  ice  sheet 
from  which  slabs  84-3  and  84-4  were  obtained.  Scale  subdivisions  on  photo¬ 
graphs  are  1  mm. 


these  crystals  extend  through  nearly  the  entire 
thickness  of  ice,  20.5  cm),  and  the  substantial  in¬ 
crease  in  the  cross-sectional  dimensions  of  crystals 
with  increasing  ice  thickness.  Such  structures, 
both  vertical  and  horizontal,  are  indistinguishable 
from  those  observed  in  normal  congelation  sea  ice 
(Gow  and  Weeks  1977,  Weeks  and  Ackley  1982). 
The  interpenetrated,  feathery  appearance  of  crys¬ 
tals  in  horizontal  sections  is  also  a  characteristic 
feature  of  congelation  sea  ice.  Though  difficult  to 
see  at  the  reduced  magnification  of  the  photo¬ 


graphs  presented  in  Figure  4,  layers  of  brine  inclu¬ 
sions  are  clearly  evident  in  the  thin  sections  them¬ 
selves.  The  bulk  salinity  of  this  ice  was  5.3%o. 

At  the  time  this  ice  was  being  sampled  for  struc¬ 
tural  and  salinity  measurements,  two  slabs,  each 
measuring  46x46  x20.8  cm  thick,  were  removed 
from  the  sheet  and  transferred  to  a  -30 °C  cold- 
room  in  preparation  for  dielectric  measurements 
(Arcone  and  McGrew  1985).  These  measurements 
were  obtained  to  determine  the  effect  on  the  di¬ 
electric  properties  of  changes  in  temperature,  sa- 
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Figure  5.  Vertical  and  horizontal  thin  section  structure  in  thermally  modi¬ 
fied  ice  slab  84-3.  Scale  subdivisions  are  I  mm. 


Unity  and  crystal  structure  of  the  slabs  as  they  are 
progressively  warmed  from  -30°  to  -2°C.  The  in- 
situ  salinity  of  the  slabs  (here  designated  as  84-3 
and  84-4)  averaged  5.3%o.  However,  brine  drain¬ 
age,  following  removal  of  the  slabs  from  the  pool, 
had  reduced  their  bulk  salinities  to  3.8%0  immedi¬ 
ately  prior  to  testing.  By  the  time  the  temperature 
of  each  slab  had  increased  to  -2°C,  the  salinities 
of  both  had  further  decreased  to  between  3.3  and 
3.2%,.  This  represents  a  real  loss  of  brine  from  the 
ice,  most  of  which  occurs  near  the  melting  point. 

The  effect  of  this  brine  loss  is  to  increase  porosi¬ 
ty,  which  in  the  case  of  slab  84-3,  amounted  to  a 


4%  increase,  as  determined  by  density  measure¬ 
ments  made  prior  to  and  at  the  end  of  testing. 
However,  the  volume  of  brine  in  the  ice  slab  actu¬ 
ally  increases  during  warm-up  because  of  eutectic 
melting  of  the  ice  around  the  brine  pockets.  This 
melting  in  turn  leads  to  progressive  enlargement  of 
the  pockets,  which  begin  to  coalesce  at  tempera¬ 
tures  of  around  -5°C\  ultimately  forming  drain¬ 
age  tubules  near  the  melting  point.  This  mobiliza¬ 
tion  and  increase  in  volume  of  brine,  together  with 
partial  drainage  of  brine  near  the  melting  point  (as 
described  above),  is  entirely  consistent  with  the  ex¬ 
planations  given  by  Arcone  and  McCirew  (Part  11, 
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I  inure  6.  Magnified  sections  o  f  ice  structure  (left)  before  and  (right)  after  thermal  modification  of  slab  S4-s 
Scale  subdivisions  are  I  mm. 


this  report)  to  account  lor  the  observed  changes  in 
the  dielectric  behavior  of  slabs  84-3  and  84-4  as 
they  underwent  warming  from  30  C  to  2  C. 

The  structural  condition  of  the  thermally  modi¬ 
fied  ice  in  slab  84-3  is  shown  in  figure  5.  I  he 
structure  departs  appreciably  from  that  observed 
in  the  sample  of  in-sit u  ice  (big.  4).  While  the  elon¬ 
gated  texture  of  crystals  in  the  thermally  modified 
ice  is  still  retained  in  vertical  sections,  crystal 
boundaries  in  the  horizontal  sections  have  become 
rounded  and  tlte  brine  inclusion  pnUcin  can  he 
seen  to  have  undergone  major  changes  when 


viewed  at  higher  magnification  alongside  sections 
from  corresponding  depths  in  the  in -situ  unmodi¬ 
fied  ice  (I  ig.6).  Crystals  in  the  unmodified  ice  ex¬ 
hibit  sharply  defined  ice  plate  brine  layer  sub¬ 
structure.  Note  that  the  brine  layers  in  this  ice  ac¬ 
tually  consist  of  strings  of  inclusions,  whereas  the 
brine  layers  in  the  thermally  conditioned  ice  have 
been  substantially  eliminated  and  replaced  by 
large  bubble-like  inclusions  and  drainage  chan¬ 
nels.  I  his  situation  applies  equally  well  to  ice  ob¬ 
tained  near  the  top  and  the  bottom  of  the  ice 
sheet.  Although  these  chances  in  civ-aal  Mihsttuc 


Figure  7.  Vertical  and  horizontal  thin  sections  of  seed¬ 
ed  ice  structure  in  slab  85-1  prior  to  microwave  testing 
and  thermal  modification.  Smallest  subdivisions  are  I 
mm. 


ture  were  induced  by  progressive  warming  of  the 
sample  in  the  laboratory,  essentially  identical 
changes  in  brine  inclusion  geometry,  salinity  char¬ 
acteristics  and  crystal  boundary  shape  were  ob¬ 
served  in  the  ice  sheet  in  the  outdoor  pool  after  it 
had  been  subjected  to  several  days  of  above-freez¬ 
ing  air  temperatures  toward  the  end  of  the  experi¬ 
ment. 

1984-85  experiments 

Experimental  measurements  of  dielectric  prop¬ 
erties  of  ice  slabs  were  extended  to  samples  taken 
from  two  of  three  ice  sheets  grown  during  the  1984- 
85  winter.  A  total  of  four  slabs  were  tested  under 
essentially  the  same  conditions  as  the  1983-84 
slabs;  that  is,  samples  were  stored  for  some  time  at 
-30°C  before  beginning  microwave  transmission 
measurements.  In  addition,  tests  were  repeated  in 
reverse  to  evaluate  the  nature  of  the  dielectric  re¬ 
sponse  as  the  slabs  were  cooled  from  -2°  to 
-30 °C.  Results  of  these  measurements  substantial¬ 
ly  confirm  the  1983-84  data  and  are  discussed  in 
detail  by  Arcone  and  McGrew  in  Part  II  of  this  re¬ 
port. 

Slabs  designated  85-1  and  85-2  were  obtained 
from  seeded  and  unseeded  parts,  respectively,  of 


Figure  8.  Vertical  and  horizontal  thin  sections  showing  original 
structure  of  unseeded  ice  m  slab  85-2.  Smutted  ccutc  cut^i/ix isum\  lire  / 
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Figure  13.  Salinity  profiles  at  differ¬ 
ent  stages  of  growth  of  ice  sheet  from 
which  slabs  85-3  and  85-4  were  har¬ 
vested. 

This  ice  sheet  was  regularly  monitored  for  salinity 
and  structure  until  19  March  1985.  In  Figure  13 
are  profiles  measured  several  days  before  and  just 
after  harvesting  sample  85-3  on  22  January  and 
just  prior  to  removing  sample  85-4  on  20  Febru¬ 
ary.  Immediately  before  testing,  the  bulk  salinity 


of  85-3  measured  5.4“/oo  and  that  of  85-4  measured 
2.4°/oo.  By  the  end  of  testing,  salinities  had  changed 
to  4.2°/oo  and  2. 3%,  respectively.  This  change 
amounted  to  a  20%  loss  of  brine  from  85-3  but 
represented  no  significant  drainage  of  brine  what¬ 
soever  from  85-4.  Whereas  85-3  was  taken  from  a 
relatively  thin  (12-cm)  brine-rich  ice  sheet  (6.2u/«, 
salinity  on  23  January)  the  ice  sheet  as  it  existed  on 
20  February  (now  24-25  cm  thick)  had  begun  to 
desalinate  from  the  effects  of  rain  that  fell  on  12 
and  13  February.  This  resulted  in  some  pooling  of 
water  on  top  of  the  ice  but  the  water  was  pumped 
off  before  it  could  freeze.  A  light  snowlall  subse¬ 
quently  led  to  upward  migration  of  brine  that 
caused  enhanced  salinities  in  the  top  2  cm  of  the 
ice  sheet,  as  demonstrated  in  the  19  February  sa¬ 
linity  profile.  However,  enhanced  salinity  at  the 
bottom  of  the  ice  sheet  (evidence  of  recent  ice 
growth),  together  with  salt  enrichment  in  the  top  2 
cm,  was  enough  to  raise  the  bulk  salinity  to  3.3"4„. 
By  the  time  85-4  was  readied  for  laboratory  testing 
in  July,  its  salinity  had  decreased  to  2.4”/,,,  mainly 
because  of  unavoidable  drainage  of  brine  during 
harvesting  of  the  sample. 

Apart  from  the  slightly  coarser  grained  nature 
of  its  crystals,  the  ice  in  85-3  (Fig.  14)  is  structural- 
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Figure  14.  Vertical  and  horizontal  thin  sections  of  structure  of  ice  in  s/ah 
85-3  prior  to  testing  and  thermal  modification  Smallest  scale  i/immm  are  I 


figure  16  Magnified  sections  of  ice  from 
slab  Hf-.i  Scute  subdivisions  are  I  mm 
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Figure  /  7.  Vertical  and  horizontal  thin  sections  of  structure  of  ice 
in  slab  85-4  prior  to  testing  and  thermal  modification.  Smaller  scale 
subdivisions  are  I  mm.  Note:  horizontal  line  at  14  cm  in  vertical  section  is 
not  a  structural  discontinuity;  it  marks  the  joining  of  two  photographs. 
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1%  identical  to  that  ot  85-1  Structure  in  the  t her 
mally  modilied  ice  (Mg.  1*)  is  also  very  similar  to 
that  tound  in  8 5  I  following  thermal  treatment,  e\ 
cept  that  the  original  substructure  ol  crsstals  in 
XVI  is  somewhat  better  preserved  than  in  KV  I . 
which  also  contained  much  larger  residual  mclti 
sums  Enlarged  structure  sections  (f  ig  16)  indicate 
the  extent  ot  brine  pocket  modification  following 
testing  and  thermal  treatment  ol  slab  1 

Ice  m  slab  X*  4  shows  verv  little  difference  in 
structure  before  or  alter  thermal  modituation  (see 
fig  P  and  IK)  I  his  is  in  accord  with  the  natural 


ly  degraded  state  of  the  ice  sheet  prior  to  sampling 
XV4  on  20  february  In  contrast  to  that  ol  KVJ 
the  dielectric  response  ol  XV4  was  also  verv  much 
suppressed,  a  situation  entirelv  compatible  with  its 
low  salimtv  characteristics 

I  xplanations  presented  bv  Arcone  and  Mctirew 
(I’art  II.  this  report)  to  account  tor  the  observed 
change  in  the  dielectric  behavior  ol  the  I4K4  H* 
saline  ice  slabs  as  thev  undergo  warming  Irom 
ft)  (to  2  (  are  entirelv  consistent  with  the  pu 
lure  ol  temperature  ion!:  oiled  brine  structure 
changes  reported  here 
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Figure  18.  vertical  and  horizontal  thin  section  of  structure  of  ice  in  slab 
85-4  after  testing  and  thermal  modification.  Scale  subdivisions  are  I  mm. 
Fine-grained,  striated  structure  at  base  of  vertical  section  and  needlehke  inclu¬ 
sions  near  middle  of  section  are  artifacts  of  thin-seciioning  process. 
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(  ONC  l.l  SIONS  AMI  RECOMMENDATIONS 

Microwave  transmission  measurements  clearly 
show  that  the  dielectric  behavior  of  saline  ice  is 
dominated  by  thermally  activated  changes  in  the 
brine  component  and  that  changes  in  relative  per¬ 
mittivity,  for  example,  are  closely  related  to  brine 
volume  changes  but  are  essentially  independent  of 
the  ice  component.  This  situation  arises  by  virtue 
ot  the  Iree/ing  characteristics  of  salt  water  that 
produce  ice  composed  of  columnar  crystals  in 
which  the  brine  is  concentrated  in  vertical  layers 
(or  sheets  ot  inclusions)  sandwiched  between 


plates  of  essentially  pure  ice.  The  activity  of  the 
brine  is  linked  closely  to  its  very  sensitive  response 
to  changes  in  temperature,  with  even  small  chang¬ 
es  resulting  in  significant  changes  in  the  size  and 
shape  of  the  brine  inclusions.  These  changes  in 
turn  control  the  mobility  and  pattern  of  migration 
of  brine  through  the  ice  matrix,  the  crystallinity  of 
which  remains  essentially  unchanged.  The  verti¬ 
cally  elongated  texture  of  congelation  crystals  in 
saline  ice  (as  in  natural  sea  ice)  ensures  that  brine 
movement  is  also  dominated  by  gravity -driven  mi¬ 
gration  in  the  vertical  dimension. 


The  need  for  detailed  measurements  of  brine  in-  ments,  would  allow  more  precise  determination  of 

elusion  geometry  and  distribution  cannot  be  over-  porosity  at  all  stages  of  thermal  modification, 

stressed.  This  has  been  partly  addressed  in  meas-  Similarly,  additional  thin  sectioning  in  both  the 

urements  made  here  on  thin  sections  of  saline  ice  vertical  and  horizontal  dimensions  would  allow 

slabs  from  the  CRREL  pool.  These  measurements  for  more  accurate  assessment  of  changes  in  the 

indicate  that  the  shapes  of  brine  pockets  are  sub-  sizes  and  shapes  of  brine  pockets  as  they  respond 

stantially  ellipsoidal,  but  that  their  sizes  can  vary  to  temperature  changes  in  the  ice.  To  assess  the  di- 

greatly,  even  in  the  different  layers  of  brine  within  electric  behavior  of  new  or  young  sea  ice  it  would 
a  single  crystal.  be  very  useful  to  harvest  slabs  under  conditions 

Any  further  evaluation  of  dielectric  properties  cold  enough  to  minimize  brine  drainage  and  en- 

by  techniques  described  here  should  include  more  sure  bulk  salinities  of  6  to  7%o  at  the  outset  of  test- 

frequent  measurements  of  density  which,  in  con-  ing. 

junction  with  salinity  and  temperature  measure- 


Part  II — Microwave  Properties 


INTRODUCTION 

The  dielectric  properties  of  sea  ice  determine  its 
emissivity  and  hence  strongly  influence  the  bright¬ 
ness  temperature.  A  program  has  been  instituted 
to  measure  microwave  emission  and  radar  back- 
scatter  of  artificially  grown  saline  ice  sheets.  It  has 
been  CRREL’s  responsibility  to  measure  the  di¬ 
electric  properties  of  these  sheets  at  several  appro¬ 
priate  frequencies  and  to  document  the  physical 
properties  and  ice  structure  as  a  function  of  the 
growth  and  decay  of  the  ice. 

Dielectric  properties  of  sea  and  saline  ice  sam¬ 
ples  and  natural  sheets  have  been  measured  at  dis¬ 
crete  frequencies  between  0.1  and  40  GHz,  but  the 
literature  is  not  extensive.  Hallikainen  (1980)  has 
made  measurements  at  600  and  900  MHz  as  a 
function  of  temperature,  as  has  Addison  (1969, 
1970)  at  100  MHz  and  less.  Bogorodsky  (1979) 
and  Morey  et  al.  (1984)  have  used  pulse  methods 
and  give  frequency  ranges  in  which  their  results 
might  apply.  Vant  et  al.  (1978)  have  measured  die¬ 
lectric  properties  as  a  function  of  frequency  for  a 
variety  of  conditions.  All  of  this  work  shows  that 
below  about  2°C ,  at  brine  volumes  of  less  than 
about  loo"/*,,  the  real  part  (or  ‘'dielectric  con¬ 
stant”)  <r  of  the  complex  relative  permittivity  <* 
generally  varies  from  between  4.0  and  5.0  near 
100  MHz  to  between  3.0  and  3.5  above  1  GHz, 
while  the  imaginary  part  (or  "dielectric  loss")  <  ’  is 
dominated  by  conductive  or  dispersive  processes 
below  10  GHz.  The  change  in  i,  with  frequency  is 
most  likely  ascribable  to  a  Maxwell-Wagner-Sil- 
lars  type  dispersion  (Sillars  1937)  as  discussed  by 
Addison  (1970)  In  none  ot  the  above-referenced 


work  was  the  structural  condition  of  the  ice  pre¬ 
sented. 

Our  objective  was  to  measure,  in  the  labora¬ 
tory,  r?  of  slab  ice  samples  at  frequencies  appro¬ 
priate  to  the  concurrent  radiometric  and  backscat- 
ter  experiments,  and  to  evaluate  these  data  in 
terms  of  the  salinity  and  structural  characteristics 
of  the  ice  slabs  as  they  were  progressively  warmed 
from  about  -32  °C  to  -2°C.  Several  sheets  were 
grown  each  year  and  seven  different  slabs  were  i  - 
moved  to  obtain  measurements  under  a  variety  of 
physical  conditions.  The  growth  history  in  terms 
of  thickness,  temperature  and  salinity  was  record¬ 
ed  for  each  ice  sheet  (Part  1).  Once  removed,  the 
slabs  were  immediately  brought  into  a  coldroom 
to  minimize  brine  drainage.  A  comparison  then  of 
the  phase  and  amplitude  of  waves  transmitted 
through  the  slabs  gave  the  forward  transmission 
coefficient  from  which  <*  could  be  calculated.  Ad¬ 
ditionally,  in-situ  measurements  of  transmission 
attenuation  using  antennas  above  and  beneath  the 
ice  sheets  were  made  to  obtain  data  for  very 
young,  brine-iich  ice.  The  data  were  then  com¬ 
pared  with  models  based  on  data  from  ice  thin  sec¬ 
tions  and  a  dielectric  mixing  theory.  This  work  is  a 
continuation  of  results  previously  presented  by 
Arcone  and  Gow  (1985)  and  Gow  (1985). 


MEASUREMENT  TECHNIQUES 

A  microwave  bridge  was  used  in  a  coldroom  tor 
these  measurements  (Fig.  19a).  The  coldroom 
temperature  was  generally  held  to  within  ■  0,1  (. 


a.  Schematic  of  microwave  bridge  with  sample  in  place. 
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Angular  Displacement  m  the  Electric  and  Magnetic  Planes 

b.  Distributions  of  amplitude  and  phase  in  the  principal  antenna  planes  across  a 
planar  interface  at  1 1  cm  from  the  transmit  horn. 

Figure  19.  Schematic  and  field  patterns  of  slab  transmission  measurement  apparatus. 

of  the  reported  measurement  temperature  during  square,  stood  vertically  during  testing  and  rested 

testing  for  any  one  slab.  Measurements  were  made  on  a  sliding  cart  to  facilitate  removal  for  an  air 

approximately  every  three  days  during  both  warm-  reference  measurement.  The  slabs  were  stored 

ing  and  cooling  cycles  to  allow  the  slabs  to  reach  horizontally  to  permit  proper  brine  drainage  when 

equilibrium.  One  arm  of  the  bridge  contained  a  it  occurred  during  warming, 

reference  section  of  coaxial  waveguide  while  the  The  distance  between  transmit  and  receive  an¬ 
other  arm  contained  two  opposing  radiators,  tennas  was  determined  by  the  need  to  utilize  plane 

which  were  either  4.8-  x  2.2-cm  open-end  wave-  wave  formulations  for  transmission  coefficients  to 

guides  (4.75  and  4.80  GHz)  or  8.6-  x  6.4-cm  stan-  interpret  t *  from  phase  and  amplitude  measure- 

dard  gain  horns  (9.50  GHz)  between  which  the  ments,  the  available  power  in  the  system,  and  the 

slabs  were  placed.  The  signals  from  both  arms  maximum  size  of  the  slabs  that  could  be  removed 

were  then  stepped  down  in  frequency  and  fed  into  from  the  ice  sheet.  A  separation  of  46  cm  placed 

a  network  analyzer  (Hewlett  Packard  Model  both  receive  antennas,  the  open-end  waveguide 

8410A)  from  which  the  relative  phase  and  gain  of  and  standard-gain  horn,  in  the  plane  wave  zone 

the  test  arm  could  be  read.  The  reading  accuracies  (i.e.  the  ‘‘far  field”  where  the  plane  wave  criteria 

were  1  "of  phase  and  0.1  dB  of  gain.  Each  radiator  of  uniform  phase  and  amplitude  over  the  receive 

could  rotate  on  its  axis  to  allow  polarization  aperture  is  satisfied)  of  the  transmit  antennas, 

changes,  and  a  microwave  absorber  surrounded  This  separation  allowed  capture  of  the  10-dB 

all  support  structures  to  reduce  unwanted  reflec-  bcamwidths  by  the  46-cm-squarc  slabs  to  mini- 


tions.  The  slabs  themselves,  approximately  46  cm  mi/e  edge  ellects,  and  permitted  a  suitable  signal 


k,  -  k0n* 


level  for  even  the  most  lossy  situation  (at  the 
warmest  temperatures  for  the  thickest  slab  at  24 
cm),  considering  the  low  power  ( <  100  mW)  avail¬ 
able  in  our  system. 

Of  equal  importance  to  the  use  of  plane  wave 
formulation  was  for  the  slabs  themselves  to  be  in 
the  plane  wave  zone  of  the  antennas.  The  standard 
mathematical  criterion  for  the  far  field  is  2 D2/\, 
(D  =  maximum  antenna  aperture  dimension,  K  is 
free  space  wavelength)  which  is  about  7  and  46  cm 
for  the  4.8-  and  9.5-GHz  systems  respectively. 
This  criterion  is  not  satisfied  by  the  9.5-GHz  sys¬ 
tem  where  a  slab  thickness  of  24  cm  would  place  a 
slab  face  at  only  1 1  cm  from  either  antenna.  How¬ 
ever,  a  physical  optics  calculation,  whereby  the 
antenna  aperture  was  divided  into  an  array  of 
point  dipole  sources  subject  to  the  appropriate 
horn  phase  and  rectangular  waveguide  TE,0  mode 
field  strength  distributions  (e.g.  Collin  1960),  re¬ 
vealed  that  3-dB  beamwidths  of  about  54°  and  32° 
had  already  formed  by  1 1  cm  in  both  the  magnetic 
and  electric  field  planes,  respectively.  At  this  dis¬ 
tance,  within  the  angles  of  ±6°  that  extend  from 
the  axis  of  the  transmit  antenna  to  the  widest  di¬ 
mensions  of  the  receiver,  both  antennas  perform 
as  point  sources  producing  a  nearly  uniform  distri¬ 
bution  of  phase  and  amplitude  over  a  slab  face 
(Fig.  19b),  thus  justifiying  the  use  of  the  plane 
wave  formulation  of  the  forward  transmission  co¬ 
efficient  subject  to  an  amplitude  correction  (dis¬ 
cussed  below)  to  account  for  refraction  collima- 
tion  of  the  slightly  diverging  rays  within  the  slab. 

Phase  and  amplitude  were  recorded  as  a  func¬ 
tion  of  temperature  and  then  compared  to  give  the 
forward  transmission  coefficient  Su,  the  complex 
ratio  of  measured  electric  field  strength  with  the 
slab  to  that  measured  without  the  slab  (i.e.  the  air 
reference).  This  quantity  was  then  used  in  the 
equation 

S,,  =  W*W,{\ -R2)/{\ ~R‘W,)  (1) 

to  solve  for  er*  such  that 


by  iteration  on  a  computer.  In  eq  1 
B'0  =  exp  (ik0d) 


n*  =  complex  index  of  refraction  =  n  '  - 
in" 

/=  v^T 

R  =  (n*~  l)/(n*+  1) 
d  =  slab  thickness. 

The  data  from  six  different  polarization  posi¬ 
tions  were  averaged  for  each  temperature,  as  we 
did  not  expect  any  anisotropy  due  to  the  random 
orientation  of  dendritic  planes  seen  between  crys¬ 
tals.  The  air  references  for  phase  and  amplitude 
were  established  for  each  polarization  at  each  tem¬ 
perature.  The  slab  faces  were  kept  normal  to  the 
transmission  axis  to  within  1  °,  as  was  the  plane  of 
rotation  of  the  radiators. 

The  quantity  (j  was  determined  mainly  by 
phase,  for  which  the  reading  accuracy  of  ±  1  ° 
gave  a  value  that  was  good  to  ±0.01.  The  phase 
readings  are  extremely  stable,  theoretically,  for 
the  slab  thicknesses  used,  varying  no  more  than  1  ° 
for  0.0  <  (’  <  0.5.  The  quantity  «r"  depended 
strongly  on  attenuation  for  which  the  reading  ac¬ 
curacy  of  ±0.1  dB  gave  a  value  that  was  generally 
good  to  ±0.003.  These  sensitivities  vary  slightly 
over  a  360°  cycle  of  phase  and  depend  also  on  slab 
thickness.  The  largest  measured  amplitude  de¬ 
crease  from  the  air  reference  for  any  one  slab  was 
-15.1  dB,  at  an  absolute  lev  1  well  above  the  noise 
range  of  the  detector. 

The  large  ratios  of  d/%,  allowed  a  multiplicity  of 
solutions,  generally  separated  by  0.6  to  1.0  for  er'. 
We  chose  those  that  tracked  nearest  3.2  at  temper¬ 
atures  below  -26  °C,  when  brine  volume  was  mini¬ 
mal  ( <  7Vo<>)  so  that  fr’  was  expected  to  be  nearly 
that  of  freshwater  ice.  The  separate  contributions 
of  dipolar  processes  and  bulk  conductivity  to  er' 
could  not  be  distinguished  by  only  phase  and  am¬ 
plitude  measurements  at  one  or  two  frequencies. 

A  gain  was  generally  measured  at  the  lowest 
temperatures  where  attenuation  from  propagation 
loss  should  be  minimal  because  of  the  minimal 
amount  of  brine  in  the  slabs.  This  gain  was  caused 
by  the  collimating  action  (dielectric  focusing)  of 
the  ice  whereby  the  divergence  of  the  radiation  is 
decreased  when  refracting  through  the  slabs.  Con- 


W,  =  exp  (ik,d) 


ko  =  2ir/\o 


sequently,  corrections  based  on  a  simplified,  geo¬ 
metric  optics  analysis  of  0.73  dB  (0.73-0.87  for 
the  thinnest  slab)  to  2.50  dB  (2.26-2.50  for  the 


Figure  20.  Microwave  attenuation  instrument  for  measur¬ 
ing  m-situ  attenuation.  The  slotted  tine  measures  reflection  loss 
of  the  submerged  transmit  antenna. 


thickest)  were  subtracted  from  all  attenuation 
readings.  This  was  a  significant  correction  at  tem¬ 
peratures  below  -6°C  where  average  attenuation 
readings  of  less  than  -7  dB  were  generally  record¬ 
ed. 

The  geometric  optics  calculation  assumed  the 
transmitting  antennas  to  be  point  sources  from 
which  rays  diverged  and  were  then  captured  by  the 
receiving  antennas,  with  and  without  the  slabs 
present.  The  gain  was  calculated  from  the  areal 
concentration  of  these  rays  after  refraction 
through  the  slabs.  For  example,  the  theoretical 
gain  caused  by  a  20.8-cm  slab  is  1.9  dB  at  er'  =3.1 
and  2.2  dB  at  er'  =  4.0.  A  15.8-cm  slab  will  give 
1.4  to  1.7  dB  of  gain  over  the  same  range  of  er'. 
These  calculations  apply  to  a  ray  bundle  confined 
to  less  than  7°  divergence  from  the  waveguide 
axis,  in  keeping  with  the  point  source  approxima¬ 
tion.  Since  each  collimating  correction  depended 
on  t / ,  the  solutions  had  to  be  recalculated  once  an 
approximate  value  of  ej  was  computed  from  the 
raw  data.  That  a  gain  was  measured  at  the  lowest 
temperatures  at  values  commensurate  with  the 
geometric  optics  prediction  implies  that  the  actual 
volume  of  ice  studied  was  extremely  small;  proba¬ 
bly  no  larger  than  20  cm1  in  cross  section  for  the 
thinnest  slabs  that  were  at  the  greatest  distance 
from  the  antennas. 

The  minimum  slab  thickness  removed  was  8.5 
cm,  as  the  intrinsic  strength  of  the  saline  ice  pre¬ 
cluded  the  taking  of  thinner  samples.  Conse¬ 
quently,  a  simple  microwave  transmission  tester 
was  also  constructed  (Fig.  20)  for  making  in-situ 
attenuation  readings.  A  slot  was  cut  in  the  ice  and 
the  transmit  antenna,  either  an  open-ended  wave¬ 
guide  (4. 5-5.0  GHz),  or  a  horn  (9-10  GHz),  was 
butted  against  the  ice  bottom.  The  receive  antenna 
was  then  located  above  the  ice  at  a  height  that  de¬ 
pended  on  the  ice  thickness.  The  butting  of  the  im¬ 
mersed  antenna  violated  the  plane  wave  assump- 


for  the  incident  wave  that  was  needed  to  interpret 
the  data  using  eq  1.  Therefore,  we  relied  on  the 
change  of  signal  amplitude  relative  to  an  air  path, 
with  corrections  made  for  the  vltage  standing 
wave  ratio  (VSWR)  of  the  submerged  horn,  and  an 
estimated  relfection  loss  at  the  air/ice  interface 
(which  was  minimal)  to  calculate  an  effective  at¬ 
tenuation  rate. 


RESULTS 

Table  1  lists  information  pertinent  to  each  slab. 
Most  of  the  slabs  were  obtained  from  seeded 
sheets.  Slabs  84-3  and  4  were  removed  at  the  same 
time,  about  24  hours  after  the  start  of  a  strong 
warming  period  in  the  weather.  Slabs  85-1  and  2 
were  also  taken  at  the  same  time.  All  slabs  were  re¬ 
moved  during  growth  and  prior  to  any  surface 
melting  except  85-4,  which  was  cut  from  the  ice 
sheet  after  it  had  experienced  one  short  period  of 
flushing  of  fresh  water  from  snowmelt  runoff. 
More  complete  documentation  of  growth  and 


Table  1.  Test  slab  data. 


Final 


Slab 

Thickness 

(cm) 

In  situ  S 

r/„) 

S  before 
testing 

S  after 
testing 

density 

(g/cm’t 

Seeded? 

84-2 

15.3 

5.0 

4.4 

4.0 

0.883 

no 

84-3 

20.8 

5.3 

3.8 

3.3 

0.884* 

yes 

84-4 

20.8 

5.3 

3.8 

3.2 

0.886 

yes 

85-1 

8.5 

8.4 

4.8 

3.1 

0.865 

yes 

85-2 

8.5 

8.4 

5.5 

4.4 

0.885 

no 

85-3 

11.8 

6.2 

5.4 

4.2 

0.866 

yes 

85-4 

24.0 

3.0 

2.4 

2.3 

0.903 

yes 

As  measured  in  neighboring  samples. 


structural  properties  of  the  several  ice  sheets  is 
given  in  Part  I. 

All  slabs  lost  brine  when  they  were  removed 
from  the  ice  sheet,  and  they  continued  to  drain  in 
the  coldroom  until  their  temperature  T  dropped 
below  -24  °C.  Brine  drainage  resumed  at  between 
-6°  and  -5  °C  during  testing.  Each  slab  was  melted 
at  the  end  of  testing  and  all  the  brine  that  had 
drained  from  it  during  testing  was  added  to  the 
meltwater  to  determine  the  salinity  S  of  the  slabs 
just  before  testing.  The  in-situ  values  of  S  in  Table 
1  are  those  for  other  slabs  removed  and  immedi¬ 
ately  sampled  for  salinity  profile  measurements  at 
the  same  time  the  slabs  for  dielectric  measure¬ 
ments  were  being  taken  from  the  ice  sheet.  Calcu¬ 
lations  of  brine  volume  uB  were  made  using  the 
equations  of  Frankenstein  and  Garner  (1967)  or 
the  tables  of  Assur  (1960)  based  on  T  and  S.  A  lin¬ 
ear  extrapolation  was  made  between  the  initial  and 
Final  values  of  S  to  calculate  after  brine  drain¬ 
age  began  during  warming.  The  values  of  uB  are 
given  in  the  figures  later  in  parts  per  thousand 
(%«)  which  may  be  converted  to  volumetric  brine 
content  by  multiplying  the  values  by  0.001 . 


1983-84  experiments  (4.75  GHz) 

The  average  measured  values  of  er'  and  e"  for 
1983-84  are  plotted  as  a  function  of  temperature 
in  Figure  21.  The  vertical  lines  indicate  the  meas¬ 
ured  range  as  polarization  was  rotated.  The  values 
within  each  range  did  not  vary  systematically  for 
each  temperature  as  polarization  was  changed. 
Slab  and  antenna  alignment  were  maintained  for 
each  polarization  measurement.  The  slab  faces 
were  flat  to  within  1  mm  although  etched  on  the 
bottom  side.  However,  the  actual  portion  of  the 
slabs  being  measured  is  estimated  at  20  cm1  or  less 
and  the  accuracy  of  laterally  repositioning  the 
slabs  between  each  polarization/air  reference  was 
no  better  than  ±2.5  cm.  Therefore,  many  differ¬ 
ent  crystal  groups  could  have  been  sampled  during 
one  polarization  scan  and  this  can  well  account  for 
the  variations  seen. 

The  largest  range  in  Figure  21  for  fr'  is  only 
±3%  for  slab  84-3  at  -2.5  °C.  The  lowest  values 
of  fr"  occur  below  -22  °C  and  are  associated  with 
the  freezing  of  almost  all  brine.  The  «r'  values  for 
slabs  84-3  and  -4  are  similar  and  slightly  lower 
than  those  of  slab  84-2  because  they  both  lost 
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more  brine  (see  Table  1)  in  going  from  the  pond  to 
the  coldroom  and  during  warming.  Both  e'  and  er" 
show  little  change  until  about  -8  °C  where  uB  starts 
to  increase  more  rapidly  with  warming  because  of 
eutectic  melting  of  the  ice  around  the  brine  poc¬ 
kets.  This  melting  in  turn  leads  to  progressive  en¬ 
largement  of  the  pockets  which  begin  to  coalesce 
at  temperatures  of  around  -5°C,  ultimately  form¬ 
ing  drainage  tubules  near  the  melting  point.  Mod¬ 
eling  studies,  presented  later,  will  show  that  this 
mobilization,  increase  in  volume  of  brine  and  par¬ 
tial  drainage  of  brine  near  the  melting  point  are 
consistent  with  the  observed  dielectric  behavior. 

Figure  22  replots  the  components  of  fr'  as  a 
function  of  eB.  The  quantity  t,  shows  a  linear  de¬ 
pendence  until  about  a  uB  value  of  50°/oo,  after 
which  the  values  do  not  rise  as  steeply  but  show 
greater  variance.  Bogorodsky  (1979)  has  presented 
a  very  similar  curve  for  100  MHz  except  that  the 
dielectric  values  are  progressively  larger  as  uB  in¬ 
creases.  The  quantity  er"  behaves  similarly  but  with 
a  greater  discrepancy  between  slabs.  Although 
slab  84-2  consistently  has  the  higher  values  of  er' ,  it 
does  not  have  consistently  higher  values  of  t"  as 
might  have  been  expected  from  a  sample  with 
higher  salinity.  Nor  are  the  e"  values  for  slabs  84-3 
and  -4  as  similar  to  each  other  as  are  their  er'  val¬ 
ues.  However,  at  the  highest  temperatures,  slab 
84-2  retained  the  most  brine  (Table  1)  and  pro¬ 
duced  the  greatest  attenuation  rate  (Fig.  23).  The 
transmission  attenuation  0  (dB/cm)  shown  in  Fig¬ 
ure  23  was  calculated  using  the  formula 

0  =  8.68(2irn"/\„). 

1984-85  experiments  (4.80  and  9.50  GHz) 

Experimental  measurements  of  dielectric  prop¬ 
erties  of  ice  slabs  were  extended  to  samples  taken 
from  two  of  three  ice  sheets  grown  during  the 
1984-85  winter.  Four  slabs  (85-1  to  4)  were  tested 
under  essentially  the  same  conditions  as  the  1983- 
84  slabs;  that  is,  slabs  were  removed  from  the  ice 
sheets  and  stored  for  some  time  at  -30  °C  before 
beginning  microwave  transmission  measurements, 
which  were  made  at  regular  intervals  during  pro¬ 
gressive  warming  of  the  slabs  from  -32°  to  -2°C. 
In  addition,  tests  were  repeated  in  reverse  to  evalu¬ 
ate  the  nature  of  the  dielectric  response  as  the 
slabs  were  cooled  from  -2°C  to  -24  °C. 
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Figure  22.  Relative  complex  permittivity  vs  calculated  brine 
volume  for  the  1983-84  slabs. 
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Figure  23.  Transmission  attenuation  vs  calculated  brine  volume 
for  the  1983-84  slabs. 


The  quantity  <*  is  plotted  vs  T  for  the  four  slabs 
in  Figures  24  and  25  at  4.80  and  9.50  GHz,  respec¬ 
tively.  The  filled  circles  are  the  initial  warming  cy¬ 
cle  and  the  open  circles  are  the  cooling  cycle.  Both 
figures  show  lower  values  of  during  the  cooling 
cycle,  as  expected  from  the  replacement  of  brine 
with  air.  The  youngest  slabs,  85-1  (seeded)  and 


85-2  (unseeded),  give  similar  t'r  values,  suggesting 
that  crystal  size  is  not  an  important  factor.  The 
oldest,  thickest  and  most  desalinated  slab,  85-4, 
shows  the  least  separation  between  warming  and 
cooling  values  of  <r'  as  it  had  the  least  loss  of 
brine.  In  addition  the  response  of  to  warming  in 
slab  85-4  was  very  much  suppressed  with  respect 
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Figure  24.  Relative  complex  permittivity  vs  temperature  for  the  1984-85  slabs,  4.80  GHz. 


to  all  other  1985  slabs,  a  situation  compatible  with 
its  low  salinity  characteristics.  Some  measure¬ 
ments,  mostly  for  slab  85-3,  show  the  cooling  cy¬ 
cle  values  below  -16°C  to  be  consistently  below 
3.18,  an  value  previously  reported  for  pure  ice 
at  -20°C  (Glen  and  Paren  1975).  These  less-than- 
pure  ice  values  are  due  to  the  presence  of  air  poc¬ 
kets  derived  from  drainage  of  some  of  the  brine 
pockets  still  remaining  in  the  original  ice  at  these 
temperatures.  The  range  of  most  measurements 
never  exceeds  about  ±  5°7o. 

Generally  slabs  85-1,  2  and  3  show  the  cooling 
cycle  values  of  «r"  to  be  lower  than  the  warming  cy¬ 


cle  values,  as  expected  from  the  loss  of  brine,  but 
the  differentiation  is  not  great.  Although  the  r,’ 
temperature  dependencies  compare  favorably  be¬ 
tween  1983-84  and  1984-85,  a  similar  comparison 
of  the  behavior  of  e"  is  distinctly  different,  either 
between  any  two  slabs  or  on  the  average.  The  val¬ 
ues  never  increase  in  the  1985  data  to  the  range  of 
0.2  to  0.4  as  they  did  in  1984.  In  fact,  slabs  85-1 
and  3  (4.80  GHz)  and  85-1  and  2  (9.50  GHz)  show 
some  very  low  values  and  erratic  behavior  in  the 
-6  to  -2°C  range.  Possible  causes  for  the  1983-84/ 
1984-85  differences  are  discussed  later. 
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Figure  25.  Relative  complex  permittivity  vs  temperature  for  the  1984-85  slabs,  9.50  GHz. 


1984-85  in-situ  experiments 

Table  2  summarizes  the  in-situ  attenuation 
measurements.  The  average  salinity  was  that  of  a 
full  thickness  sample  and  the  average  temperature 
is  the  mean  of  the  surface  and  bottom  tempera¬ 
tures.  The  loss  was  measured  relative  to  an  air 
path  whose  own  loss  left  only  about  a  50-  to  60-dB 
range  for  measurement.  A  correction  was  made 
for  reflection  losses  at  the  ice/submerged  transmit 
antenna  interface  (measured)  and  at  the  air/ice  in¬ 
terface  (estimated).  The  seven  measurements  in 
Table  2  not  in  the  noise  range  of  the  detector  vary 
between  4.1  and  8.7  dB/cm.  The  highest  loss  rate 
occurred  in  the  thinnest  ice,  which  was  also  rela¬ 
tively  warm  and  thus  had  the  highest  value  of  rH. 
These  calculated  values  of  i  H  are  generally  high 
and  complement  the  range  of  <  H  values  in  figure 


23.  By  comparison  with  Figure  23,  Table  2  shows 
that  the  attenuation  rate  ,i  will  start  to  increase 
dramatically  above  100“/IH,. 

The  9-  to  10-GH/  attenuation  rate  is  always 
larger  than  the  4.5-  to  4.8-CiHz  rate  over  the  first 
eight  values  of  Table  2.  This  contrasts  with  the 
lower  brine  volume  data  of  Figures  24  and  25, 
which  when  averaged,  will  be  shown  below  to  give 
more  comparable  i,"  values  for  the  two  frequencies 
in  1984-85.  Similar,  high  in-situ  attenuation  rates 
that  increase  with  frequency  for  young  ice  have 
also  been  reported  by  Bogorodsky  (1979).  Attenu¬ 
ation  rales  due  only  to  d.c.  conduction  currents 
are  very  similar  at  both  frequencies  (discussed 
later)  and  so  cannot  be  responsible  for  this  attenu¬ 
ation  increase  with  increasing  frequency.  The  di¬ 
electric  mixing  model  discussed  in  the  next  section 
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Table  2.  In  situ  loss  measurements,  1985. 


Transmission 


Sheet 

number 

Thickness 

(cm) 

Frequency 

(GHz) 

: 

!| 

Loss 

(dB/cm) 

T 

ro 

S 

ru)  _ 

Esl  i  j 
f'U) 

2 

7.1 

4.5 

41.6 

5.9 

-4.5 

9.1 

lOl 

2 

7.1 

10.0 

43.5 

6.1 

-4.5 

9.1 

101 

2 

9.1 

4.8 

37.5 

4.1 

-3.0 

9  2 

149 

2 

9.1 

9.0 

>52.6 

>5.8 

-3  0 

9.2 

149 

3 

5.8 

4.8 

29.6 

5.1 

-2.2 

8  4 

183 

3 

5.8 

9.0 

50.5 

8.7 

-2.2 

8.4 

183 

3 

6  9 

4.8 

36.9 

5.3 

-2. Of 

6.7 

161 

3 

6.9 

9.15 

51.8 

7.5 

-2.0 

6.7 

161 

3 

12.0 

4.8 

>52.2 

>4.4 

-2.9 

5.8 

97 

3 

12.0 

9.0 

>52.2 

>4.4 

-2.9 

5.8 

97 

3 

24.0 

4.8 

>52.2 

>2.2 

-l.l 

3.2 

146 

3 

24.0 

8.9 

>52.2 

>2.2 

-1.1 

3.2 

146 

*  Brine  volume  calculated  from  the  equations  of  Frankenstein  and  Garner  (1967)  and  the  tables 
of  Assur  (I960)  using  average  values  of  temperature  and  salinity, 
t  Surface  temperature  is  estimated. 


does  predict  an  attenuation  increase  with  fre¬ 
quency,  but  requires  a  bulk  d.c.  conductivity  term 
to  give  the  right  attenuation  order  of  magnitude. 

DISCUSSION 
Comparison  of  data 

Figure  26  plots  the  averages  of  based  on  all 
values  occurring  within  increments  of  5%,,.  Aver¬ 
ages  were  taken  only  where  at  least  three  data 
points  in  an  increment  occurred.  The  plots  show 


that,  during  both  warming  and  cooling  cycles,  r 
was  slightly  greater  at  4.80  GHz  than  at  9.50, 
which  is  consistent  with  past  dispersional  models 
of  sea  ice  (Addison  1970,  Vant  et  al.  1978)  and  the 
one  discussed  later.  Figure  26  also  shows  other 
values  reported  at  4  and  10  GHz  (Bogorodsky  and 
Khoklov  1975,  Vant  et  al.  1978;  all  values  given  in 
Vant  et  al.  1978)  for  samples  for  which  brine  vol¬ 
umes  were  calculated  from  the  given  values  of  T 
and  S.  There  is  reasonable  agreement  up  to  about 
60%o  at  which  point  many  of  our  data  points  seem 
to  merge.  This  is  misleading,  however,  as  our 


f  igure  26.  >,  i  s  i  H  f'W  comparisons  heheecn  ituia  averages  lor  /V,VT  XJ 
and  1984  <95  slabs,  and  with  data  <>!  previous  research  tor  < . 


measurements  above  507*,  at  9.50  GHz  varied  be¬ 
tween  3.6  and  4.3,  a  variation  similar  to  that  of 
Bogorodsky  and  Khokiov  (1975). 

Figure  26  also  shows  that  at  any  one  value  of  rB, 
the  1984-85  cooling  cycle  value  of may  be  2  to 
8%  lower  than  the  warming  cycle  value.  A  conver¬ 
sion  of  the  "before”  and  "after"  salinities  (Table 
1)  into  brine  volumes  (at  -2°C)  gives  a  maximum 
percentage  gain  in  air  volume  of  about  3  to  4% 
due  to  the  drained  brine.  For  example,  the  salinity 
changes  in  slab  85-3  translate  to  a  density  change 
during  testing  from  0.906  to  0.872  g/cm’  (actual 
tinal  value  =  0.865  g/cm')  or  about  3%  gain  in  air 
volume.  The  average  observed  decrease  in  at, 
for  example,  <■„  =  40°/x  is  4.0%  for  all  1985  slabs 
and  4.3%  for  slabs  85-1,  85-2  and  85-3,  which 
drained  the  most  brine. 

This  close  correlation  between  decrease  in 
and  gain  in  air  volume  during  warming  suggests 
that  modifications  in  the  brine  structure  from 
thermal  effects  do  not  significantly  affect  if  This 
is  consistent  with  the  thin-section  studies,  which 
show  retention  of  individual  brine  inclusions  in 


the  thermally  modified  ice,  although  without 
strong  alignment  along  dendritic  planes.  This  cor¬ 
relation  is  also  consistent  with  the  comparison  ot 
1984  averages  with  those  for  the  1985  cooling  cy¬ 
cles  (Fig.  26)  which  is  close  because  slabs  84-3  and 
-4  had  values  5-7%  lower  than  the  1985  slabs’ 
warming  values.  These  lower  values  of  >,  could  be 
partially  caused  by  some  loss  of  brine  structure 
(i.e.  discrete  inclusions)  due  to  warm  weather,  but 
more  probably  to  a  loss  in  brine  volume.  The  air 
temperature  at  the  time  of  removal  of  these  slabs 
was  near  7°C  (from  a  previous  overnight  low  of 
-2.9°C)  and  much  brine  was  lost  upon  extraction 
and  transport  of  the  slabs  to  the  coldroom.  This 
loss  of  brine  was  estimated  using  the  salinity  data 
of  Table  1  and  used  to  estimate  the  initial  densities 
of  the  slabs  before  dielectric  testing.  The  values 
obtained  were  utilized  for  the  modeling  discussed 
in  the  next  section. 

Figure  27  compares  <"  data  averages  from 
1983-84  and  1984-85  with  ranges  of  values  from 
Bogorodsky  and  Khokiov  (1975)  and  Vant  et  al. 
(1978).  The  1984-85  data  are  only  for  the  warming 
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Figure  27.  t ,  vs  i  fl  (“/<*,)  comparisons  between  data  averages  of  1983-84  and 
1984-85  slabs  and  with  data  of  previous  research  for  if  Open  circles  are  maxi¬ 
mum  values  measured  for  indicated  slabs  at  4.80  GHz. 


Table  3.  Saline  ice  structural  parameters. 


Brine  mi  tuuons 

1  v  t'rugt'  Ih'fHh  nt  /  hntlnlt  ' 
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C  b  tLUiil 
rat  to  fS  l)  ) 
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b 
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(mm)  (S  O  ) 
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anat\zed 

(i  mi 
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: 
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rHI  RM  MOD 

85-2 

0  460  (0  pi) 

1  46  0  05 

0  40  0  03 

0  91 

(0  41) 
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12 

0  W  id  lit 
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No 

0  442  (0  199) 
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M  4  l  |M  1  Ni 

THFRM  MOD 

84-3,4 

0  417  (0  188) 

0  35  0.04 

0  19  0  0) 

0  9* 
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4  20  (1 
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11  H)  ill] 

84-3,4 
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0  439  (0  197) 

0  74  0  09 

0  21  0.01 

0  VI 

(0  15)e 
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4  20  i> 
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a — Range  pertains  only  to  inclusions  measured  in  a  random  sampling  Actual  limits  slightls  greater  lor  less  tor  minimal 
b— 20  crystals  surveyed 

c— Drainage  structure  had  developed;  could  not  trace  dendrites 
d— C  an  change  with  depth 

e— Mainly  small  circles  along  vertical  lines,  some  very  large  pits  (3  A  mm) 
f— Dendrites  could  not  be  distinguished 


cycle  to  facilitate  the  comparison.  The  plots  show 
the  distinction  between  the  1983  84  and  1984-85 
slabs  in  t  B  dependency  above  about  20°/,x,.  The 
1983-84  data  and  a  few  isolated  maximum  values 
measured  in  1984-85  are  consistent  with  Vant  et 
al.  (1978)  af4  GHz.  The  9.50-GHz  values  at  t  B  < 
50%,,  are  slightly  lower  than  those  at  4.80  GHz  for 
the  1985  slabs;  this  is  not  consistent  with  the  theor¬ 
etical  model  discussed  below,  which  predicts  slight¬ 
ly  higher  values  at  9.50GHz  due  to  the  brine  whose 
loss  peaks  around  9  GHz.  However,  the  order  ot 
magnitude  is  consistent  with  the  model.  The  agree¬ 
ment  with  Vant  et  al.  (1978)  at  10  GHz  is  poor,  al¬ 
though  there  is  some  overlap  There  is  no  overlap 
with  the  data  of  Bogorodsky. 

Possible  explanations  for  the  difference  be¬ 
tween  the  1983-84  and  1984-85  values  of above 
brine  volumes  of  about  20%o  were  investigated  but 
are  inconclusive.  The  difference  seems  related  to 
bulk  conductivity,  which  is  discussed  later.  The 
phase  and  gain  meters  of  the  network  analyzer 
were  all  found  to  be  in  calibration  after  the  1985 
testing,  although  we  cannot  discount  temporary 
malfunctioning.  (The  extreme  dependence  of  >, 
upon  phase  precluded  any  possibility  of  an  uncali¬ 
brated  gain  meter  affecting  the  validity  of  the 
calculations.)  Structural  thin  sections  (Table  3)  ex¬ 
amined  prior  to  testing  revealed  the  1983  84  brine 


pockets  to  be  vertically  longer  and  narrowet  than 
those  seen  in  the  1984  85  slabs,  and  this  mav  he 
related  to  greater  conductivity.  However,  these 
sections  were  taken,  not  from  the  slabs,  bin  Irom 
adjacent  samples  because  we  did  not  want  to  do 
turb  the  slabs  prior  to  testing.  All  thin  sections  ol 
thermally  modified  ice  were  taken  from  the  'labs 
themselves.  It  would  have  been  desirable  to  re 
move  additional  slabs  for  simultaneous  warming 
and  thin  sectioning  to  examine  structural  changes 
occurring  during  testing,  but  this  was  not  done  in 
order  to  minimize  disturbing  the  ice  sheet  lor 
other  experiments 

Some  general  dit  lei  cnees  between  t  he  198  3  84 
and  1984-85  slabs  are  evident  (  ompared  with 
1984,  density  is  lower  in  slabs  8^  |  and  3  (both  ot 
which  had  very  low  values  ol  <,  neat  4  t  |,  age 
(or  thickness)  is  less  m  slabs  85  1.  2  and  3,  and  the 
temperature  ot  initial  freezing  lot  all  sheets  w,i' 
much  higher  in  1988  (freezing  was  initialed  at 
about  20  (  in  1984  hut  only  at  about  10  (  in 
1985)  I  he  modeling  data  (below  I  show  that  the  • 
values  below  ■  „  no  ,.  can  be  m.ikhed  lor  the 
1984  85  slabs  by  assuming  a  brine  volume  dcpcit 
deni  bulk  conductivity  that  never  exceeds  0  03 
S  m  (siemens  meter),  whereas  values  up  to  0  0b 
S  m  are  generally  required  loi  1481  ,84  slabs 
Iheretore.  a  greater  contumitv  ol  sondiisiive 
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channel'  (i.e.  pen  eabilii > )  ma>  have  existed  in 
ihe  1484  84  slabs,  possible  due  lo  their  older  con 
dmon. 

Analytical  modeling 

•\n>  theor>  of  dielectric  properties  in  the  \  HI 
through  microwave  range  should  account  tor  ll  a 
decrease  in  with  increasing  frequency.  2)  an  in 
crease  in  both  <,  and  >,  with  increasing  1)  at 
least  two  loss  mechanisms,  with  total  loss  (.,  )  in 
creasing  as  frequency  decreases  into  the  HI  range, 
and  4)  vers  high  microwave  losses  at  high  values 
ot  i  H.  with  attenuation  increasing  with  (requeue v 
Ihe  dielectric  mixing  model  given  in  \  ant  et  al 
(1478)  and  discussed  in  the  appendix  meets  these 
criteria  with  the  following  exceptions  It  as  tre 
quenev  decreases  toward  and  below  100  Mil/.  ■, 
eventually  increases  an  impossibility  lor  the 
known.  Unite  values  ol  sea  ice  condustivitv  (  I  li\s 
sen  et  al  |4'4,  I  unco  14^4).  and  2)  the  predicted 
losses  at  high  values  ol  H  are  not  a'  high  as  oh 
served  Ihe  error  lies  in  the  lack  ot  an  additional 
loss  mechanism,  namely  bulk  sonduciiv itv  u, 
which  reflects  the  permeabilitv  ot  sea  ice  In  addi 
non,  our  thin  section  studies  reveal  a  much  cl 1 1 let 
ent  brine  inclusion  geornetrv  than  was  assumed 
and  derived  bv  \  ant  et  al  |I4',K)  Ml  dieleetiic 
mixing  models  consider  on l v  unconnected  indu 
sums 

Ihe  \  ant  et  al  ( 1 4"’8 >  model  assumes  that  the 
large  and  increasing  values  ot  that  occur  as  tie 
queries  decreases  toward  100  Mil/  result  Irom  an 
assumed  angular  displacement  "  ot  1'  4^  ot  the 

long  axes  ot  ellipsoidal  brine  pockets  Mom  the 
(settled)  directum  ot  wave  propagation  I  hi' 
angle  "  would  also  be  the  displacement  ot  the  ice 
ctvstal  c  axes  Itom  horizontal  Ihe  ellipsoid'  ate 
assumed  to  have  circular  cto"  sections  ti  e  are 
prolate  spheroids),  the  diameters  ot  which  are 
much  shoitei  than  that  ot  the  long  axis  so  ih.it  a 
moderate  value  ol  "  is  needed  to  permit  stioitg 
coupling  ol  the  hoti/outal  electric  Meld  into  the 
brine  ellipsoids  Although  tins  model  will  give 
higher  loss  near  UK)  MHz.  ev  ent  nail  v  • ,  will  de¬ 
crease  as  Irequetlcv  decreases  below  KMi  MHz,  at 
anv  "  value  I  Ins  decrease  is  not  compatible  wnh 
the  behavior  ot  a  conductive  material 

Prior  to  moditvmg  tilts  model,  measiitement' 
were  made  on  enlargements  ot  the  thin  section 
photographs  shown  m  Part  I  Brine  pocket  geom 
etrv  m  both  the  horizontal  and  vertical  planes  was 
examined  tor  several  slabs  and  the  data  are  sum 
marized  m  I  able  1  Below  about  belli  in  the  an 
seeded  ice  and  about  I  cm  m  the  ceded  he.  tin 


dendritic  planes  Iwhich  arc  normal  n>  :  tie  ,  axe-i 
seen  in  vertical  ttim  sections  -tn*w  iuiIv  marginal 
1  ill  w  hie  tl  as  cr  ag’e  I  '  <  u  U-s-  'mii,  vc!’i,.il  I  Ik 

hon/onlal  I  Inn  section  -how  Mu  mo-. 
the  brine  inclusions  to  r aiige  Horn  oiciil.n  to  chip 
Inal  (long  /»  axes  parallel  to  the  planes,  ,  .m-  pet 
pendic ulat )  In  the  hotizontal  plane  then,  the 
hi  me  pockets  ate  genetallv  ellipsoidal  a-  well  .t' 
being  ellipsoidal  m  the  vertical  plane  m  which  the 
v  et  t  leal  J  axis  is  -still  r  tie  longest  I  her  el  or  e .  our  et 
hpsoids  contain  three  unequal  axial  length-  with 

tl  I  ho  ellipsoidal  georneliv  in  the  hotizonial 
plane  gives  the  same  qualitative  result-  tot  ■  a 
\  ant  el  al  tl4"X|  obtained  and  .tlso  predict-  that 
■  must  eventuallv  decrease  with  decreasing  (re¬ 
queue  \ .  -it  -in  v  c  a  I  tie  ol  "  I  tiei  et  oi  e  .i  -oiiductiv 
tc  loss  let  in  is  still  required 

I  tie  \  ant  el  al  i  |4~xi  model  o  theietoie  modi 
lied  tn  selling  "  0  and  adding  a  bulk  d  ,  ,oii 

ductivitv  loss  term  (  i  .  ,  -  wtiere  i-  tre 

queue v  m  rad  s and  -  »  x'  •  In  I  mi  t  on 

duel iv it v  values  ranging  Mom  o  tH  to  o  IKH  s  m 
lia  v  e  been  reported  lor  sea  h  e  (  I  hi  seen  et  ai  I  4' 4 
1  unco  i4"’4t,  values  sutMcierit  to  produce  the  Iocs 
es  observed  hv  \  ant  et  al  1 1 4"’8 1  at  IBM  MHz  >hu|k 
has  been  measured  mdirectlv  through  mierpreta 
noils  ot  Schlumberger  sounding  curves  tor  ne 
sheets  grown  in  I48h  (Areone  I4X'i  Ihe  values 
ranged  from  0*1  (t  200  . )  to  (>  1  i  H  2  .  t 
tor  treshly  growing  ice  wuh  a  great  deal  ot  .aria 
bihtv 

figure  28  compares  the  4  xt)  (.Hz  e  per itnental 
data  with  the  revised  model  leq  \4  Appendix  Al. 
emplovmg  the  data  <  ’  t  I  able  P  depolat  izat  ion  so 
elltcienis  lor  ellipsoid'  ot  htee  unequal  axial 
length'  tt  is  (so  i  n  144*'),  and  bt  me  -  undue  t  iv  it  v  and 
dtelecttic  values  Mom  Stogtvn  and  I  )e'.it  cent 
II4XM  lire  hgures  list  the  mean  axial  lengths  p.i 
each  slab,  the  t.tiios  ot  wludi  wne  used  to  ,om 
pule  the  depolat  iz.il  ion  c  oel  !  i.  lent  -  t  fi.i!  ate  -tbo 
eiceu  Ihe  ratio'  >t  import. u.d  .u.  >  ,.  .md  . 
which  were  computed  Mom  the  average  lu  me  m 
elusion  duileiisioii-  I  lu  .  j . . .  i  r  1 1 1  r  v  ,  •  o  i  v  or,  u. 

I  .d'li  1  is  itii  mi-, in  ol  !  tic  individual  ratio-,  o'  ’tn 
horizontal  .tv  -  .md  not  'lu  ’.uio  .■■  'In  un.ii 

lengths 

I  he  models  give  good  .plan t  it -U  i v  i  ,oiiip,ii  soii> 
with  'ot  slat’s  '4  t  .md  4  nut  toi  ;  lie  .oolitic 
cVcIes  of  tat'--  x'  2  and  1  I  hi'se  food  ,1'inp.iM 
soils  aie  toi  slabs  it  u  had  .iheudv  desaluiaied  ap 
pti'ciablv  .tint .  [hellion,  piot'.it'lv  ti.nt  weaker 
eleclrh.il  inter. h  lion-  I’eiweeu  inclusion-  lu  .ill 
c.i'i's  both  ttieoiellc.il  .md  txpi’  Mu  nt.it  wanning 

V  C  ll  d.ll  .1  O  *  t  cl  .  I  !  .  ,  ;  ;  .  I  '  ’  1  ,  :  '  -  A  :  '  ’  I’  I  ’  I  1  i 
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28.  Comparisons  of  experimental  data  with  predictions  of  dielectric  mixing  theory  presented  in  the  Appendix.  Brine  di- 
parameters  taken  from  Stogryn  and  Desargenl  (1985);  brine  pocket  depolarization  coefficients  taken  from  the  curves  of  Osborn 
brine  inclusion  model  ellipse  is  in  the  horizontal  plane  with  b  the  direction  of  the  dendritic  planes;  depolarization  coefficient  P/a  = 
*  P. );  a,  b,  c  are  semi-axis  lengths;  ahr  =  brine  conductivity. 
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f  igure  29.  Model  results  for  < vs  brine  volume  at  4.80  GHz. 


-8°  and  above  -24  °C  in  the  rate  of  decrease  of  fr' 
with  decreasing  temperature.  The  change  near 
-8°C  results  from  changes  in  brine  electrical  prop¬ 
erties  and  the  change  at  -24  °C  is  associated  with 
the  precipitation  of  sodium  chloride  dihydride 
(Assur  1960). 

Slab  85-3  (Fig.  28)  shows  greater  theoretical  dis¬ 
crepancy  between  warming  and  cooling  cycles 
than  does  85-2,  reflecting  the  slightly  greater  loss 
in  brine  suffered  by  slab  85-3  and  lower  final  den¬ 
sity.  Only  small  differences  between  warming  and 
cooling  cycle  depolarization  coefficients  were  cal¬ 
culated  for  each  slab,  eliminating  this  as  a  likely 
explanation  for  the  differences.  The  similarity  of 
these  axial  ratios  and  of  the  initial  salinity  of  both 
slabs  before  warming  produce  very  similar  warm¬ 
ing  curves  for  i, . 

Improved  quantitative  comparisons  with  <r"  are 
achieved  by  adding  a  bulk  d.c.  conductivity  term, 
except  in  the  case  of  slab  85-3.  Linear  and/or  Win- 
sauer  type  exponential  (Ward  and  Fraser  1967) 
empirical  dependencies  of  ohu|k  upon  rB  were  as¬ 
sumed  in  one  case  as  was  done  by  Thyssen  et  al. 
(1974).  These  empirical  conductivity  terms  do  not 
explicitly  state  any  dependency  upon  permeability; 
the  dependency  is  implicit  in  the  brine  volume  pro¬ 
portionality.  The  permeability  factor  is  the  most 
likely  cause  of  differences  in  >j'  between  the 
1983-84  and  the  1984-85  slabs.  The  lack  of  consis¬ 
tent  differentiation  between  warming  and  cooling 
cycle  experimental  values  in  Figures  28b  and  28c 
implies  that  brine  volume  is  not  as  strong  an  influ¬ 
ence  upon  dielectric  contributions  to  if  as  it  is 
upon  fr' . 

Figure  29  replots  the  4.80-GHz  modeling  results 
for as  a  function  of  r„.  The  cooling  curves  show 


linear  dependency  to  55%o,  the  highest  brine  vol¬ 
ume  attained  in  this  cycle.  The  warming  curves  are 
all  linear  to  between  40  and  60%o  and  then  show  a 
slight  decrease  in  slope  before  resuming  another 
linear  increase.  The  modeling  results  for  if  in  this 
brine  volume  range  will  naturally  be  linear  (or 
nearly  so)  because  of  the  assumed  linear  or  non¬ 
linear  dependence  of  conductivity  upon  rH. 

The  Vant  et  al.  (1978)  model  and  our  revision 
both  predict  increasing  loss  with  increasing  micro- 
wave  frequency,  but  do  not  give  the  large  values 
we  and  others  have  observed  at  rB  >  lOO0/™  if  i;hulk 
is  not  included.  This  low  loss  is  demonstrated  in 
Figure  30  which  gives  theoretical  r*  values  and  the 
associated  attenuation  rates  ,1  for  a  hypothetical, 
high  brine  volume  sea  ice  model  (based  on  the 
brine  inclusion  parameters  of  slab  85-3)  at  the  two 
microwave  frequencies  studied.  The  values  of  o’ 
for  "bulk  =  0  arc  much  lower  than  observed  (Table 
2).  However,  if  0.2  <  ohu|k  <  0.5  S/m  at  brine  vol¬ 
umes  >  100"/™,,  then  d  (Fig.  31)  would  fall  in  the 
range  of  values  observed  in  Table  2  while  still 
maintaining  the  2-  to  4-dB  difference  between  4.80 
and  9.50  GHz.  Note  that  Figure  31  shows  attenua¬ 
tion  rates  due  only  to  d.c.  conduction  to  be  nearly 
equal  at  both  frequencies,  so  that  d.c.  conduction 
cannot  be  responsible  for  an  increase  of  attenua¬ 
tion  with  frequency.  The  high  values  of  ,>  observed 
in  Table  2  therefore  suggest  that  bulk  conductivity 
does  have  an  exponential  dependency  on  brine 
volume  that  increasingly  dominates  attenuation  as 
<  )t  increases  above  100" . 

Although  the  structural  and  salinity  characteris¬ 
tics  of  the  saline  sheets  simulate  natural  condi¬ 
tions,  removal  of  slabs  or  any  other  type  of  sample 
is  always  accompanied  b>  brine  drainage,  as  ev  i - 
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Figure  30.  Dependency  of  <  *  and d  upon  temperature  for  a 
hypothetical  model  of  high  hrine  volume  sea  ice.  The  data 
are  intended  to  model  the  behavior  of  freshly  growing  ice.  The  as¬ 
sumption  of  ahu,k  =  0  shows  the  inadequacy  of  the  model  in 
predicting  observed  attenuation  rates  Id). 
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Figure  31.  Attenuation  rate  vs  dielectric 
constant  for  various  values  of  conduc¬ 
tivity  a  at  4.75  and  9.50  GHz  < broken 
curve).  Only  by  o  =  I  S/m  does  the  9.50GHz 
rate  differ  noticeably  from  that  for  4.  75  GHz 
so  that  d.c.  conduction  cannot  be  responsible 
for  the  increase  in  attenuation  rates  seen  with 
increasing  microwave  frequency. 


deneed  by  the  density  values  given  in  Table  1.  Fig¬ 
ure  30  shows  that  at  any  given  brine  volume,  use 
of  the  in-situ  sea  ice  density  of  0.926  elevates  all 
modeling  fr'  values  for  warming  cycles  by  about 
3%  and  fr"  values  by  2. 5%  or  less  (compare  with 
Fig.  28c  for  slab  85-3).  Therefore,  despite  the  loss 
of  brine,  the  dielectric  permittivity  values  are  be¬ 
lieved  to  simulate  well  the  properties  as  they  exist 
in-situ. 

A  modification  of  this  model  (as  well  as  other 
scattering  models)  based  on  inclusion  interaction 
was  examined  in  hopes  of  obtaining  a  better  match 
of  fr'  to  the  warming  curves  of  slabs  85-2  and  3.  A 
measured,  average  dendrite  length  of  2.41  mm  and 
a  width  of  0,03  mm  were  used  for  the  ellipsoidal 
axial  lengths  in  the  h  and  (  directions  respectively 
to  account  for  the  close  spacing  •  *  the  inclusions 
along  dendritic  planes.  These  dimensions  give  a 
very  small  value  to  l\  (  =  0.01 )  which  results  in  < , 
giving  a  good  match  for  T  <  -9  but  with  errors 
approaching  30,ro  in  excess  of  our  data  as  7  ap¬ 
proaches  -2  .  Unfortunately.  <,  exceeds  our  data 


by  an  order  of  magnitude  and  is,  thus,  far  in  ex¬ 
cess  of  all  other  reported  data  at  these  frequencies. 
A  theory  based  on  multiple  scattering  in  a  dense 
medium  (Tsang  and  Kong  1980)  gives  even  less 
loss  than  that  predicted  by  linear  independent  scat¬ 
tering  theory  (e.g.  Kerr  1965)  which  gives  the  right 
frequency  dependence  but  values  that  are  at  least 
an  order  of  magnitude  too  low  .  More  recent  devel¬ 
opments  by  Stogryn  (1985)  may  hold  more  prom- 
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The  i,  results  for  the  seven  saline  slabs  at  two 
frequencies  all  show  a  dependence  on  cB  that  is  lin¬ 
ear  to  about  50“/™,  and  then  nonlinear  above  50"/, 
with  progressively  greater  variance.  C  omparison 
of  warming  and  cooling  values  revealed  that  the 
replacement  of  brine  with  air  could  have  lowered 
i,  by  as  much  as  6°'<>  at  any  one  salue  of  >  while 
changes  due  to  any  redistribution  ol  brine  inclu- 
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sions  were  probably  not  significant.  The  1983-84 
slab  and  1984-85  in-situ  data  for  e"  are  consistent, 
as  they  generally  show  an  increase  with  increasing 
iiB.  These  data  were  consistent  with  only  a  few  iso¬ 
lated  values  of  t "  for  the  1984-85  slabs.  Slabs  re¬ 
moved  from  seeded  and  unseeded  sections  of  the 
same  ice  sheets  (these  slabs  differ  mainly  in  the 
size  of  their  columnar  crystals)  showed  similar  di¬ 
electric  properties.  Modeling  studies  suggest  that 
in-situ  dielectric  values  at  eB  <  100%<,  were  proba¬ 
bly  no  more  than  about  3%  greater  than  the 
warming  values  measured  for  the  removed  slabs. 

Coupling  between  brine  inclusions  and  the  elec¬ 
tric  field  vector  is  given  for  causing  the  behavior 
of  ej  at  microwave  frequencies  as  has  been  pro¬ 
posed  and  formulated  by  others  (Tinga  et  al.  1973, 
Vant  et  al.  1978).  However,  this  coupling  must  re¬ 
sult  mainly  from  an  average  inclusion  elongation 
in  the  horizontal  plane  as  was  observed,  not  from 
an  angular  displacement  of  the  dendritic  planes 
relative  to  the  (usually  vertical)  direction  of  wave 
propagation  as  has  been  previously  hypothesized 
(Vant  et  al.  1978)  at  least  in  first-year  sea  ice.  The 
two  concepts  are  mathematically  similar  because 
the  coupling  effects  of  vertically  tilted  ellipses 
come  from  their  projection  onto  the  horizontal 
plane,  and  the  theory  does  not  specify  absolute  el¬ 
lipsoid  size,  only  semi-axis  ratios.  Consequently, 
there  is  much  leeway  in  the  theory  to  account  for 
brine  inclusion  shape  and  orientation,  while  ice 
and  brine  properties  are  well  documented  (Glen 
and  Paren  1975,  Stogryn  and  Desargent  1985).  No 
present  theory  considers  any  effect  due  to  the  ab¬ 
solute  vertical  length  of  the  brine  inclusions. 

The  choice  of  a  particular  geometry  for  the 
brine  pockets  necessitated  neglect  of  interactions 
between  inclusions  higher  than  first  order  in  the 
mixing  theory  (Tinga  et  al.  1973).  This  cannot  be 
correct  for  first-year  ice,  where  brine-filled  inclu¬ 
sions  located  between  the  dendritic  plates  of  the 
crystals  can  be  spaced  closer  than  an  inclusion 
length  and  often  appear  interconnected.  Conse¬ 
quently,  the  model  best  fits  the  fr'  data  in  Figure 
28  for  the  more  desalinated  or  multi-year  type  ice, 
which  contains  more  brine  inclusions  that  have 
drained.  This  does  not  suggest  that  electrical  per¬ 
meability  in  the  horizontal  plane  can  decrease  with 
thermal  modification,  only  that  there  is  less  di¬ 
polar  interaction  between  the  inclusions.  There¬ 
fore,  the  role  of  conductivity  should  be  more  thor¬ 
oughly  investigated,  especially  as  a  function  of  ag¬ 
ing  which  affects  permeability.  This  is  being  done 
on  artificial  sheets  grown  during  1986  (Arcone 
1987).  Field  measurements  of  conductive  and  mi¬ 


crowave  properties  of  anisotropic  ice  would  be  of 
immense  help  in  clarifying  the  role  of  various  loss 
mechanisms.  Campbell  and  Orange(l974),  Kovacs 
and  Morey  (1978),  and  Arcone  and  Delaney  (1979) 
discuss  VHF  observations. 

The  mixing  theory  used  in  this  research  can  also 
be  used  to  predict  some  frequency  dependence  as 
was  measured  by  Vant  et  al.  (1978).  This  theory 
predicts  that  er'  will  slowly  decrease  with  increas¬ 
ing  frequency  throughout  the  microwave  and  mil¬ 
limeter  wave  region,  which  is  consistent  with  the 
comparison  of  e'  at  9.50  and  4.80  GHz  seen  in  Fig¬ 
ure  26.  Above  10  GHz,  tj  will  generally  lie  between 
3.2  and  4.0  at  brine  volumes  up  to  100%o.  Values 
of  t"  at  frequencies  above  20  GHz  cannot  be  pre¬ 
dicted  with  this  theory  because  scattering  occurs 
when  in-situ  wavelengths  approach  the  circumfer¬ 
ential  dimension  of  a  brine  pocket,  which  can  be 
several  millimeters  judging  by  the  data  in  Table  3. 
Quasi-optical,  free-space  transmission  techniques, 
such  as  done  here,  or  Fabret-Perot  interferometry 
could  be  used  for  measurements  at  higher  frequen¬ 
cies,  but  highly  collimated  beams  would  sample  a 
smaller  volume  of  ice  than  did  our  system,  and 
slabs  would  have  to  be  much  thinner  than  8  cm  to 
obtain  significant  values  of  er".  Sensing  of  a  small 
sample  volume  is  thought  to  cause  the  ranges  of 
values  seen  about  our  data  points.  These  ranges 
may  enlarge  at  higher  frequencies  given  the  diffi¬ 
culties  of  precise  realignment  over  a  long  period  of 
storage  and  measurement. 

The  current  experiments  were  limited  to  congel¬ 
ation  ice.  However,  future  measurements  should 
include  tests  on  granular  or  simulated  frazil  ice. 
Frazil  ice  occurs  typically  in  the  top  layer  of  arctic- 
sea  ice  and  its  dielectric  properties  may  significant¬ 
ly  differ  from  those  of  congelation  ice  because  of 
major  differences  in  the  structural  makeup  of  the 
two  ice  types.  Since  frazil  ice  is  a  major  compo¬ 
nent  of  the  earliest-formed  sea  ice,  its  microwave 
dielectric  properties  are  of  obvious  importance  to 
the  remote  sensing  community. 
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APPENDIX  A:  DIELECTRIC  MIXING  MODEL  OF  SEA  ICE 

This  model  was  originally  proposed  by  Taylor  (1965,  1966),  extended  by  Tinga  et  al. 
(1973)  and  adapted  to  sea  ice  by  Vant  et  al.  (1978).  Tinga  et  al.’s  model  solves  for  the  dielec¬ 
tric  permittivity  of  a  mixture  of  two  or  three  different  dielectrics  consisting  of  a  host  medium 
(<,),  which  contains  any  density  of  inclusions  in  the  form  of  concentric  ellipsoids  (tj,  fi)  that 
are  smaller  than  any  wavelength  imposed  upon  the  medium.  Therefore,  static  theory  (La¬ 
place’s  equation)  is  used  to  find  the  resulting  electric  field  within  the  ellipsoids  (and  conse¬ 
quently,  the  relative  permittivity)  due  to  an  impressed,  constant  electric  field  across  the  ma¬ 
terial.  The  individual  volumes  of  the  ellipsoids  are  of  no  consequence  because  of  the  wave¬ 
length-size  criterion.  The  parameters  of  importance  are  the  relative  volume  of  the  inclusions, 
their  axial  ratios  as  expressed  through  the  ellipsoidal  depolarization  coefficients,  and  the  ori¬ 
entation  of  the  ellipsoid  axes  relative  to  the  impressed  electric  field. 

The  application  to  sea  ice  as  adapted  by  Vant  et  al.  (1978)  gives  the  mixture  dielectric  con- 


Yl  <■(*>  -<i)  -An 

'av  fl  V,  (-(P1/P,)P,(e1-«  ,)  +  #>,(<, -i,)  +  «,] 

where  V,  and  V \  are  the  volume  fractions  of  the  host  medium  (ice)  and  ellipsoidal  inclusions 
(brine)  respectively.  The  quantities  P,,  Pi  are  the  principal  depolarization  coefficients  for  the 
ellipsoidal  directions  parallel  to  the  applied  electric  field  and  depend  only  on  the  semi-axis 
ratios.  These  quantities  are  equal  for  any  particular  direction  because  the  sample  volume  of 
the  host  material  surrounding  the  brine  inclusion,  and  in  which  the  internal  electrical  field  is 
calculated,  is  assumed  to  be  bounded  by  an  ellipsoidal  surface  confocal  with  the  inclusion  el¬ 
lipsoid.  According  to  Tinga  et  al.  (1973),  the  use  of  this  additional  ellipsoidal  surface  ac¬ 
counts  for  inclusion  interactions  to  first  order  and  thus  permits  application  of  eq  Al  to  ma¬ 
terials  with  a  high  density  of  inclusions. 

Equation  Al  can  be  reiterated  to  include  air  (f2)  inside  the  host  medium  (<av)  to  give  a  new 
value  to  <av  because  air  bubbles  are  not  confocal  with  brine  ellipsoids.  The  calculations  used 
to  produce  the  curves  in  Figure  28  assumed  spherical  air  bubbles  as  observed  in  the  thin  sec¬ 
tions.  A  further  refinement  (not  utilized  here)  would  be  to  assume  some  of  the  air  pockets  arc 
ellipsoidal  because  of  any  drained  brine.  However,  this  would  offer  only  a  minor  change  in 
f *  values. 

Vant  et  al.  (1978)  offer  further  refinements  by  utilizing  the  extraordinary  ray  formula 
from  optical  crystallography  to  compute  the  dielectric  constant  <av.  „  for  sea  ice  in  which  the 
long  axis  of  the  brine  ellipsoid  is  tilted  at  an  angle  »  from  vertical  (i.e.  the  direction  of  inci¬ 
dent  wave  propagation): 


where  a  and  ft  refer  to  the  long  and  short  semi-axis  lengths  of  the  vertical  ellipsoid,  respec¬ 
tively  (c  is  the  second,  short  semi-axis  of  Vant  et  al.  (1978)  where  c  =  ft).  They  then  assume 
random  orientation  in  azimuth  of  the  electric  field  with  respect  to  the  ellipsoids,  to  give  an 
effective,  overall  dielectric  constant  (5ff  such  that 
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It  is  proposed  here  that  P,(  =  P:)  of  eq  A1  correspond  to  either  the  b  or  c  ellipsoid  direc¬ 
tion  in  the  horizontal  plane  and  that  eq  A3  be  supplemented  with  the  loss  term 

— '  ffbulk^ ue° 

where  ffbu)k  is  the  bulk  d.c.  conductivity  of  the  ice.  The  reader  is  referred  to  Timco  (1979)  for 
a  possible  derivation  of  abulk  from  the  constituent  material  properties.  In  addition,  eq  A2  is 
disregarded  for  the  main  bulk  of  first-year  type  ice  because  0  is  on  the  average  close  to  0°. 
Equation  A3  is  then  simply  altered  to 

feff  —  (fav,c  fav,b)^^  —  *  ^bulk^^o  (A4) 

to  account  for  the  random  orientation  in  the  horizontal  plane. 

In  keeping  with  the  ellipsoidal  notation  of  Vant  et  al.  (1978),  our  a  axis  is  in  the  vertical  di¬ 
rection  while  b  and  c  are  in  the  horizontal  directions.  The  depolarization  coefficients  Plb  and 
Plc  in  the  b  and  c  directions,  respectively,  were  taken  from  the  graphs  of  Osborn  (1945)  be¬ 
cause  closed  form  integral  solutions  are  not  available  for  cases  where  the  three  axes  are  un¬ 
equal.  Use  of  these  tables  requires  a  simple  redefinition  of  a,  b  and  c  if  the  measured  axes  do 
not  always  decrease  in  length  in  the  same  sequence.  The  depolarization  coefficients  used  in 
eq  Al  are  normalized  by  4 ir  as  shown  by  Tinga  et  al.  (1973). 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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